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Heat  Transfer  of  Supercritical  Carbon  Dioxide 
Flowing  in  a Cooled  Horizontal  Tube 


Douglas  A.  Olson 
Process  Measurements  Division 
National  Institute  of  Standards  and  Technology 
Gaithersburg,  MD  20899 


Abstract 

We  report  measurements  of  heat  transfer  coefficients  of  flowing  supercritical  carbon  dioxide 
(7.38  MPa  critical  pressure,  3 1 .1  °C  critical  temperature)  in  a cooled  horizontal  tube.  The  tube 
was  10.9  mm  ID,  was  cooled  over  274  cm,  and  had  an  adiabatic  entrance  section  of  55.9  cm. 
Cooling  was  accomplished  by  flowing  cold  water  countercurrent  to  the  carbon  dioxide  in  an 
annular  gap  between  the  inner  tube  (12.7  mm  OD)  and  an  outer  tube  (16.6  mm  ID).  This  set  a 
convective  boundary  condition  similar  to  what  would  be  encountered  in  a counterflow  heat 
exchanger.  The  four  experimental  parameters  found  to  affect  the  heat  transfer  coefficients,  and 
their  ranges  tested,  were:  operating  pressure  (7.409  MPa  to  13.0  MPa),  CO2  mass  flow  rate  (1.1 
kg/min  to  5.0  kg/min),  rate  of  heat  removal  by  cooling  (1780  W to  6220  W),  and  CO2  average 
temperature  (4.3  °C  below  to  24.2  °C  above  the  pseudocritical  temperature).  The  Reynolds 
number  range  at  the  CO2  average  temperature  was  63  250  to  291  700.  Conditions  which  resulted 
in  the  highest  heat  transfer  coefficients  were  pressure  close  to  the  critical,  high  mass  flow  rate, 
and  temperature  close  to  the  pseudocritical.  For  CO2  average  temperatures  lower  than  1 °C 
below  the  pseudocritical,  higher  rates  of  cooling  increased  the  heat  transfer  coefficient,  whereas 
for  temperatures  higher  than  1 °C  below  the  pseudocritical,  higher  rates  of  cooling  reduced  the 
heat  transfer  coefficient.  The  measured  Nusselt  numbers  were  predicted  to  within  experimental 
uncertainty  by  the  Krasnoshchekov-Protopopov  correlation  for  supercritical  flow,  when  the  CO2 
average  temperature  was  above  the  pseudocritical.  This  correlation  was  developed  for  constant 
heat  flux  heating  instead  of  cooling.  No  existing  correlation  was  found  to  predict  the  measured 
Nusselt  number  when  the  CO2  average  temperature  was  lower  than  the  pseudocritical.  The 
measured  heat  transfer  coefficients  were  higher  than  those  predicted  by  the  constant  property 
correlations  throughout  the  range  of  parameters  tested. 

Key  words:  carbon  dioxide;  cooled  supercritical  flow;  counterflow  heat  exchanger;  heat  transfer; 
horizontal  tubes;  supercritical  fluids;  turbulent  flow 


Certain  commercial  equipment,  instruments,  materials,  or  software  are  identified  in  this  paper  to  foster 
understanding.  Such  identification  does  not  imply  recommendation  or  endorsement  by  NIST,  nor  does  it  imply  that 
the  materials  or  equipment  identified  are  necessarily  the  best  available  for  the  purpose. 


1 


1. 


Introduction 


Carbon  dioxide  has  been  proposed  as  a “natural  refrigerant”  in  various  kinds  of  refrigeration,  air 
conditioning,  and  heat  pump  cycles  (Lorentzen  and  Pettersen,  1993).  For  equivalent  cycle 
efficiencies,  carbon  dioxide  offers  much  reduced  ozone  depletion  potential  and  global  warming 
potential  than  chlorofluorocarbons  and  hydrofluorocarbons.  It  is  also  non-toxic  and  non- 
flammable. However,  system  designs  will  have  to  be  adapted  to  the  characteristics  of  carbon 
dioxide,  and  the  cycle  efficiency  will  be  important  in  determining  whether  the  systems  will  be 
commercially  viable.  In  contrast  to  cycles  with  conventional  refrigerant  fluids,  the  pressure  in  a 
carbon  dioxide  cycle  will  be  above  the  critical  during  the  heat  rejection  process.  Hence  a 
supercritical  cooling  operation  would  replace  the  gas  condenser.  Very  little  is  known  about  the 
heat  transfer  of  supercritical  fluids  undergoing  cooling. 

In  this  work  we  have  measured  the  heat  transfer  performance  of  supercritical  carbon  dioxide  that 
is  cooled  while  flowing  in  a horizontal  tube.  (The  carbon  dioxide  critical  pressure,  Pc,  is  7.38 
MPa  and  the  critical  temperature,  Tc,  is  31 .1  °C).  In  this  experiment,  carbon  dioxide  at  pressures 
higher  than  the  critical  flowed  through  a 10.9  mm  ID  tube  and  was  cooled  by  water  flowing 
countercurrent  to  it  over  2.74  m of  its  length.  The  flow  of  CO2  was  turbulent  throughout  the 
tube.  The  primary  experimental  parameters  that  we  varied  were  system  pressure,  flow  rate, 
cooling  rate,  and  temperature.  We  report  how  the  heat  transfer  varies  with  the  experimental 
parameters,  and  compare  the  data  to  correlations  in  the  literature.  This  work  is  a companion 
experiment  to  the  one  we  reported  earlier  (Olson  and  Allen,  1998).  In  that  experiment,  we 
measured  heat  transfer  with  supercritical  carbon  dioxide  in  the  same  tube,  but  the  tube  was 
heated  instead  of  cooled  by  flowing  water.  Taken  together,  these  results  are  a comprehensive 
description  of  the  heat  transfer  in  a supercritical  fluid  that  exchanges  heat  with  a secondary  fluid. 

Although  heat  transfer  in  supercritical  fluids  has  been  studied  by  many  investigators,  as 
summarized  by  Hall  and  Jackson  (1978),  nearly  all  of  this  work  has  been  for  flows  undergoing 
heating.  As  pointed  out  by  Pitla  et  al.  (1998),  the  few  correlations  that  exist  for  cooling 
situations  show  significant  inconsistencies.  Supercritical  flows  are  complex  because  of  the  rapid 
variations  in  transport  and  thermophysical  properties  that  occur  near  the  critical  point  (Fig.  1 ; 
data  from  Span  and  Wagner,  1996,  and  Vesovic  et  al.,  1990).  When  heat  transfer  measurements 
are  compared  to  those  expected  for  a constant  property  fluid  with  these  fluid  properties,  the  heat 
transfer  can  be  enhanced  or  degraded,  depending  on  flow  or  heating  conditions  that  have  no 
effect  in  a constant  property  fluid.  Our  earlier  work  (Olson  and  Allen,  1998)  showed  that  a 
correlation  developed  for  supercritical  flows  undergoing  constant  wall  heating  can  be  used 
successfully  when  the  boundary  condition  is  a convective  one.  These  results  showed  that  both 
the  heat  flux  and  the  system  pressure  cause  the  heat  transfer  of  a supercritical  fluid  to  deviate 
from  that  of  a constant  property  fluid. 

2.  Description  of  experimental  apparatus 

2.1  Flow  loops 

Experiments  were  conducted  with  the  counterflow  heat  exchanger  test  section  placed  in  the  heat 
transfer  facility  shown  schematically  in  Figs.  2 and  3.  The  facility  consists  of  two  closed  fluid 
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loops,  one  for  carbon  dioxide  and  one  for  water,  which  join  at  the  test  section.  The  flow  loops 
were  filled  with  instrument  grade  carbon  dioxide  (99.99  % purity)  and  deionized  water, 
respectively. 

Flow  was  circulated  through  the  carbon  dioxide  loop  by  a centrifugal  pump  with  a magnetically 
coupled  drive,  which  had  a maximum  capacity  of  72  L/min  (at  no  pressure  drop).  A variable 
speed  drive  on  the  motor  of  the  pump  was  used  to  control  the  flow  rate.  Carbon  dioxide  was 
added  to  the  system  with  a 500  mL  syringe  pump,  which  also  controlled  the  operating  pressure. 
Flow  rate  was  measured  with  a coriolis  flow  meter.  After  exiting  the  flow  meter,  the  carbon 
dioxide  entered  an  “electric  heater”  heat  exchanger,  which  raised  its  temperature  to  the  desired 
test  section  inlet  temperature,  about  89  °C  at  the  maximum.  The  electric  heater  consisted  of  5.5 
m of  19.1  mm  OD  316  stainless  steel  tubing,  wrapped  on  the  outside  with  a 9 kW  heater  coil. 

The  heater  coil  was  powered  by  a 208V AC,  50  amp  SCR  controlled  power  supply.  Setting  the 
power  to  the  coil  controlled  the  temperature  of  the  carbon  dioxide  leaving  the  heater  section. 
Upon  exiting  the  heater  section,  the  carbon  dioxide  entered  the  test  section,  which  was  a 
counterflow  heat  exchanger  with  cold  water  flowing  in  the  annular  gap  outside  of  the  tube 
containing  the  carbon  dioxide.  After  leaving  the  test  section,  the  carbon  dioxide  was  cooled 
down  to  ambient  temperature  by  another  counterflow  heat  exchanger,  referred  to  as  the  “gas 
cooler”,  utilizing  laboratory  cold  water  as  the  cold  sink.  This  heat  exchanger  contained  5.5  m of 
12.7  mm  OD  tubing  for  the  heat  exchange  area.  The  carbon  dioxide  then  reentered  the  pump. 

The  system  volume  was  about  2 L.  The  tubing  connecting  the  test  section  outlet  to  the  gas 
cooler  was  19  mm  OD  316  stainless  steel,  and  the  tubing  from  the  heater  to  the  test  section  inlet 
connection  was  a Teflon-lined,  stainless  steel  braided  hose.  All  of  the  remaining  tubing  was  12.7 
mm  OD  316  stainless  steel.  All  components  in  the  flow  loop  were  rated  to  a pressure  of  13.8 
MPa  or  higher. 

The  water  flow  loop  was  similar  in  design  to  the  carbon  dioxide  flow  loop.  The  water  was 
circulated  with  a turbine  pump,  rated  for  19  L/min  at  40  m of  head.  The  flow  rate  was  controlled 
with  the  pump  bypass  valve.  Flow  rate  was  also  measured  with  a coriolis  flow  meter.  After 
exiting  the  flow  meter,  the  water  entered  the  “primary  cooler”.  The  primary  cooler  was  a 
counterflow  heat  exchanger  containing  5.5  m of  19  mm  OD  stainless  steel  tubing,  with 
laboratory  chilled  water  (5  °C)  as  the  cold  sink.  The  flow  rate  of  the  chilled  water  to  the  primary 
cooler  could  be  controlled  from  0 to  1 5 L/min.  The  flow  loop  water  entered  the  “secondary 
cooler”  after  leaving  the  primary  cooler.  This  heat  exchanger  also  contained  5.5  m of  19  mm 
OD  stainless  steel  tubing,  with  a 60/40  mixture  of  ethylene  glycol/water  flowing  countercurrent 
to  the  water.  The  temperature  of  the  glycol/water  mixture  was  controlled  by  pumping  it  from  a 
recirculating  cooler  (3  kW  capacity  at  0 °C,  controllable  temperature).  Setting  the  flow  rate  of 
the  chilled  water  to  the  primary  cooler,  in  combination  with  the  temperature  of  the  glycol/water 
mixture  in  the  recirculating  cooler,  controlled  the  temperature  of  the  water  entering  the  test 
section.  The  carbon  dioxide  in  the  test  section  then  heated  the  water.  An  air-operated  piston 
pump  was  used  to  fill  the  flow  loop  and  maintain  pressure  when  the  system  was  operating.  A 
back-pressure  regulator  prevented  over-pressure  of  the  system  due  to  thermal  expansion  of  the 
water.  The  water  loop  was  rated  to  pressures  up  to  2100  kPa,  and  was  operated  at  about  700  kPa. 
All  connecting  tubing  was  19.1  mm  OD,  316  stainless  steel. 

The  electric  heater  of  the  carbon  dioxide  loop  was  insulated  with  2.5  cm  of  basalt-fiber  insulation 
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overlaid  with  1.3  cm  of  fiberglass  pipe  insulation.  The  connecting  tubing  between  the  electric 
heater  and  the  test  section  was  insulated  with  1 .3  cm  of  fiberglass  pipe  insulation.  The  gas 
cooler  and  all  remaining  connecting  tubing  in  the  carbon  dioxide  loop  were  insulated  with  1.3  cm 
of  neoprene  insulation.  The  primary  and  secondary  coolers  of  the  water  flow  loop  were  insulated 
with  1 .3  cm  of  neoprene  insulation.  All  of  the  connecting  tubing  of  the  water  flow  loop  was 
insulated  with  1.3  cm  of  fiberglass  insulation. 

2.2  Test  section 

The  experimental  test  section  is  shown  in  Fig.  4.  It  was  a counterflow  heat  exchanger  consisting 
of  a 12.7  mm  OD,  10.9  mm  ID  inner  tube  mounted  concentric  to  a 19.1  mm  OD,  16.6  mm  ID 
outer  tube.  Carbon  dioxide  flowed  from  left-to-right  inside  the  inner  tube,  and  water  flowed 
from  right-to-left  in  the  annular  space  (1.9  mm  radial  gap)  between  the  two  tubes  (as  per  Fig.  4). 
The  tubes  were  mounted  at  both  ends  in  O-ring  compression  fittings  that  sealed  the  carbon 
dioxide  and  permitted  axial  movement  due  to  thermal  expansion.  These  fittings  were  coupled  to 
swaged  fittings  in  which  the  flow  loop  inlet  and  outlet  tubing  was  connected.  The  length  of  the 
inner  tube  was  348  cm.  The  first  55.9  cm  of  tubing  was  surrounded  by  insulation  only,  and 
hence  there  was  no  heat  transfer.  The  outer  tube,  from  centerline  to  centerline  of  the  connection 
fittings,  extended  over  the  next  274.3  cm  of  the  inner  tube.  The  final  17.8  cm  of  the  inner  tube 
length  also  was  insulated  only.  The  length-to-diameter  ratio  of  the  entrance  section  was  5 1 
(based  on  the  inner  tube  inner  diameter),  and  the  length-to-diameter  ratio  of  the  heated  section 
was  25 1 . 

The  inlet  and  outlet  temperatures  of  the  carbon  dioxide  and  the  water  to  the  test  section  were 
measured  in  thermowells  (made  from  commercially  available  3/4  inch  “T”  fittings)  that 
permitted  the  insertion  of  the  thermometers  and  allowed  the  fluid  to  flow  past  the  thermometers. 
The  test  section  outlet  thermowells  were  located  a minimum  of  15  cm  from  the  last  90°  bend  to 
the  test  section,  which  provided  sufficient  tube  length  to  mix  the  fluids  to  a uniform  temperature. 
Pressure  taps  were  located  in  the  carbon  dioxide  flow  loop  at  the  thermowells.  Pressure  taps 
were  located  in  the  water  flow  loop  a short  distance  upstream  of  the  test  section  inlet  and 
downstream  of  the  outlet.  The  inlet  pressure  of  both  loops  was  connected  to  absolute  pressure 
transducers,  and  the  pressure  difference  between  the  inlets  and  outlets  was  measured  with 
differential  pressure  transducers.  The  entire  test  section  was  mounted  on  a rigid  aluminum  I- 
beam.  Slotted  holes  with  Teflon  shims  provided  axial  movement  of  the  test  section  relative  to 
the  I-beam,  which  resulted  from  thermal  expansion.  The  test  section  was  insulated  with  1 .3  cm 
thick  fiberglass  pipe  insulation,  as  were  all  the  thermowells,  thermometers,  and  connecting 
tubing  from  the  flow  loops. 

2.3  Instrumentation 

We  measured  the  mass  flow  rate,  test  section  inlet  temperature,  test  section  outlet  temperature, 
test  section  inlet  pressure,  and  test  section  differential  pressure  for  both  the  carbon  dioxide  and 
water  sides  of  the  test  section.  Both  of  the  mass  flow  rates  were  measured  with  coriolis  mass 
flow  meters.  The  carbon  dioxide  flow  meter  had  a range  of  0 kg/min  to  10  kg/min,  with  a 
manufacturer’s  stated  relative  uncertainty  of  ± 0.15  % of  reading,  ± 0.68  kg/h  zero  stability.  The 
water  flow  meter  had  a range  of  0 kg/min  to  30  kg/min,  with  a manufacturer’s  stated  relative 
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uncertainty  of  ± 0.10  % of  reading,  ±0.163  kg/h  zero  stability.  (For  this  and  other 
manufacturer’s  stated  uncertainties,  we  assume  the  values  given  are  normal  distributions  and  that 
they  have  about  a 95  % “confidence  interval”.  Hence  the  “standard  uncertainty”  is  l/2  the  stated 
value.  See  Taylor  and  Kuyatt  [1994]  for  a detailed  description  of  the  NIST  uncertainty  policy). 

We  measured  fluid  temperatures  with  platinum  resistance  thermometers,  6.4  mm  diameter  and 
15  cm  long.  The  probes  were  calibrated  by  the  manufacturer  and  had  a combined  repeatability 
and  hysteresis  uncertainty  of  ± 0.05  °C.  The  probe  measuring  the  carbon  dioxide  inlet 
temperature  was  immersed  to  a depth  of  5.1  cm;  the  remaining  probes  were  immersed  to  a depth 
of  7.6  cm.  Self-heating  and  stem  conduction  errors  were  calculated  to  be  less  than  ± 0.001  °C. 
The  carbon  dioxide  pressure  at  the  test  section  inlet  was  measured  with  a quartz  crystal  pressure 
transducer  with  a 41  MPa  full  scale  range.  The  combined  repeatability,  linearity,  and  hysteresis 
uncertainty  of  the  transducer  was  ±4.1  kPa.  Water  absolute  pressure,  water  differential  pressure, 
and  carbon  dioxide  differential  pressure  were  measured  with  variable-reluctance  pressure 
transducers.  These  transducers  each  had  a combined  repeatability,  linearity,  and  hysteresis 
uncertainty  of  ± 0.25  % of  full  scale;  or  an  absolute  uncertainty  of  5200  Pa,  860  Pa,  and  86  Pa, 
respectively  for  the  water  absolute  pressure,  water  differential  pressure,  and  carbon  dioxide 
differential  pressure. 

Signals  from  the  thermometers  and  variable  reluctance  pressure  transducers  were  multiplexed 
through  an  automated  scanner  and  measured  with  a digital  voltmeter.  Signals  from  the 
flowmeters  were  measured  with  a frequency  counter.  The  signal  from  the  quartz  crystal  pressure 
transducer  was  converted  to  an  ASCII  character  string  by  the  manufacturer’s  electronics.  Relay 
switching  transients  and  voltmeter  A/D  conversion  introduced  negligible  uncertainty  in  the 
measured  quantities.  Measurement  of  the  signals  was  controlled  by  a personal  computer.  Raw 
signals  and  converted  parameters  were  stored  on  the  hard  drive.  Some  of  the  parameters  were 
displayed  on  the  video  monitor  to  assist  in  monitoring  and  operating  the  experiment.  Heat 
transfer  performance  parameters  were  calculated  at  the  completion  of  the  experiments. 

3.  Description  of  experiments  and  analysis  techniques 

3 . 1 Experiments  conducted  and  procedure 

Seven  “data  sets”  of  experiments  were  conducted  and  their  conditions  are  summarized  in  Table 
1.  Each  set  was  performed  on  a separate  day.  The  major  conditions  that  were  varied  for  the  set 
are  shown  in  bold  typeface.  The  pressure  throughout  the  test  section  was  above  the  critical  for 
all  experiments.  For  data  sets  1 through  4,  the  carbon  dioxide  pressure  was  fixed  at  a constant 
value,  while  the  carbon  dioxide  flow  rate  and  the  cooling  heat  transfer  rate  were  varied. 
Typically,  four  cooling  rates  were  tested  at  five  flow  rates  each.  For  these  four  data  sets  we  also 
held  the  average  of  the  carbon  dioxide  inlet  and  outlet  temperatures  (Tb)  approximately  constant 
at  1 to  2 °C  above  Tm.  The  pseudocritical  temperature,  Tm,  is  the  temperature  at  the  operating 
pressure  where  the  specific  heat  reaches  its  maximum  (slightly  above  the  critical  temperature), 
and  as  the  results  which  follow  will  show,  it  is  an  important  parameter  in  determining  the  heat 
transfer  performance.  For  four  of  the  79  settings,  the  Tb-Tm  temperature  difference  exceeded  2 
°C.  This  occurred  at  high  cooling  rates  and  low  CO2  flow  rates,  and  was  due  to  the  minimum 
temperature  possible  in  the  water  loop  (9  °C). 
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In  data  set  5,  the  CO2  pressure  was  varied  at  constant  flow  rate  and  several  values  of  the  cooling 
rate.  The  pressure  was  decreased  to  within  0.3  % of  the  critical.  The  average  CO2  temperature 
was  kept  approximately  constant,  which  resulted  in  changes  in  Tt~Tm  as  the  pressure  varied.  In 
data  set  6,  the  pressure  was  held  constant  at  7.8  MPa,  the  flow  rate  was  kept  constant,  and  the 
average  temperature  was  varied  for  several  cooling  flow  rates.  Finally,  in  data  set  7,  we  varied 
the  flow  rate  and  heating  rate  at  constant  pressure,  this  time  with  the  average  temperature  colder 
than  Tm. 

Once  a setting  was  achieved,  the  instruments  were  scanned  and  signals  stored  multiple  times 
(usually  11).  All  of  the  measured  and  calculated  values  for  each  scan  of  the  instruments,  each 
experimental  setting,  can  be  found  as  Excel  2000  files  on  the  diskette  in  the  back  cover  on  the 
report.  Altogether  we  achieved  142  experimental  settings  and  sampled  1564  data  points. 
Measured  parameters  are  listed  first,  followed  by  calculated  performance  parameters,  and  finally 
standard  uncertainties  of  the  most  important  parameters.  For  the  first  8 settings  of  data  set  3,  the 
water  pressure  drop  was  not  measured.  It  was  estimated  based  on  the  flow  rate  and  density. 

We  filled  the  carbon  dioxide  loop  in  a consistent  manner  to  achieve  a high  purity  of  carbon 
dioxide.  Beginning  with  air  at  atmospheric  pressure  in  the  system,  the  loop  was  pressurized  with 
carbon  dioxide  gas  at  the  saturation  pressure  of  the  supply  bottles  (about  5.7  MPa).  This  mixture 
of  air  and  carbon  dioxide  was  vented  down  to  a pressure  slightly  above  atmospheric,  preventing 
back-filling  with  air.  The  procedure  of  filling  with  carbon  dioxide  and  venting  to  atmosphere 
was  performed  a minimum  of  three  times.  Assuming  uniform  mixing  of  the  initial  charge  and 
the  added  carbon  dioxide  during  each  fill,  this  reduced  the  air  concentration  to  a volume  fraction 
of  less  than  1 x 10'6.  After  the  fourth  fill  of  carbon  dioxide  gas,  we  filled  the  syringe  pump  with 
liquid  carbon  dioxide  from  a bottle  containing  a dip  tube.  The  pump  dispensed  the  liquid  into  the 
flow  loop.  The  syringe  pump  was  re-filled  from  the  bottle  and  the  liquid  dispensed  into  the  loop 
until  the  loop  pressure  began  to  increase,  indicating  all  the  vapor  space  was  filled  with  liquid. 
Using  the  syringe  pump  in  a constant  pressure  mode  set  the  system  pressure.  Prior  to  circulating 
the  carbon  dioxide  and  starting  heat  transfer,  the  syringe  pump  was  set  at  about  half-full.  This 
allowed  it  to  both  dispense  and  withdraw  fluid  as  the  specific  volume  of  carbon  dioxide  in  the 
flow  loop  changed  due  to  temperature  changes. 

To  generate  heat  exchange  between  the  carbon  dioxide  and  the  water  in  the  heat  exchanger  test 
section,  both  circulating  pumps  were  turned  on.  The  gas  cooler  in  the  CO2  loop  was  turned  on  to 
maintain  a high  density  in  the  fluid  entering  the  pump.  Turning  on  and  adjusting  the  power  to 
the  electric  heater  increased  the  temperature  of  the  carbon  dioxide.  Heat  was  removed  from  the 
water  and  its  temperature  was  adjusted  by:  (1)  turning  on  the  lab  chilled  water  flow  to  the 
primary  cooler;  and  (2)  setting  the  temperature  of  the  recirculating  cooler  which  pumped  coolant 
to  the  secondary  cooler.  In  some  cases,  settings  which  required  high  water  temperatures  allowed 
us  to  shut  off  the  secondary  cooler.  For  some  settings  at  low  cooling  rates,  we  could  turn  off  the 
primary  cooler  and  use  the  secondary  cooler  by  itself.  The  four  major  independent  parameters, 
which  allowed  setting  the  conditions  for  an  experiment,  were  the  power  to  the  CO2  electric 
heater,  the  recirculating  cooler  temperature,  the  carbon  dioxide  flow  rate,  and  the  carbon  dioxide 
pressure.  We  adjusted  the  heater  power  and  the  recirculating  cooler  temperature  until  the  carbon 
dioxide  average  temperature  and  heat  transferred  achieved  the  desired  values. 
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Achieving  a new  steady  setting  after  any  of  the  four  parameters  were  changed  required  about  15 
min.  Large  changes  in  temperature  required  more  time,  and  small  temperature  changes  or  small 
flow  rate  changes  required  less  time.  We  defined  this  stability  as  a condition  when  the  changes 
of  the  experimental  parameters  or  measured  quantities  with  respect  to  time  were  small  enough  so 
that  negligible  errors  were  introduced  into  the  calculated  performance  parameters  (<  0.1%). 

Flow  and  pressure  steady  state  occurred  within  seconds  of  changing  the  set  point,  whereas 
thermal  steady  state  required  temperatures  of  the  electric  heater,  recirculating  cooler,  and  test 
section  to  stabilize.  Observing  the  temperature  changes  at  the  fluid  thermometers  indicated 
thermal  stability.  A drift  rate  of  0.1  °C/min  was  low  enough  to  introduce  negligible  error  in  the 
performance  parameters. 

3.2  Analysis  to  determine  heat  transfer  coefficient 

We  characterized  the  heat  transfer  performance  of  the  supercritical  carbon  dioxide  by  calculating 
the  heat  transfer  coefficient  from  the  experimental  data.  The  heat  exchanger  test  section  had 
three  major  modes  of  heat  transfer:  convection  from  the  hot  carbon  dioxide  to  the  inner  wall; 
radial  conduction  through  the  tube  wall;  and  convection  from  the  inner  tube  to  the  water.  The 
basic  principle  of  the  data  analysis  was  to  calculate  the  overall  heat  transfer  performance  and 
then  subtract  out  the  tube  conduction  and  water  convection  to  get  the  carbon  dioxide 
performance.  Axial  conduction  in  the  tube  and  two  fluids  was  insignificant.  Conduction  out  the 
insulation  to  the  ambient  room  was  a few  percent  and  was  accounted  for  in  calibration. 

Heat  transfer  in  a counterflow  heat  exchanger  can  be  expressed  by  the  “log  mean  temperature 
difference”  equation  which  defines  the  overall  heat  transfer  coefficient,  U (Rohsenow  and  Choi, 
1961): 


Q — UAht&TiM , 

where 


(1) 


A Tlm  = 


ATl-AT0 


In 


A7^ 

ATn 


Here,  Q = 
Aht  = 

U = 
T = 
ATlm  - 
atl  = 


A To  = 


total  heat  transferred  between  the  water  and  carbon  dioxide  over  the  entire 
heat  exchanger  length; 

heat  transfer  area  between  the  water  and  carbon  dioxide,  based  on  the  inner 

tube  inside  radius; 

overall  heat  transfer  coefficient; 

temperature,  either  water  (H2O)  or  carbon  dioxide  (CO2); 
log  mean  temperature  difference; 

Tco2  - Th20  at  the  “Z”  location  of  the  heat  exchanger,  which  is  the  CO2  outlet 
and  the  water  inlet;  and 

Tco2  - Th20  at  the  “0”  location  of  the  heat  exchanger,  which  is  the  CO2  inlet 
and  the  water  outlet. 
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If  the  overall  heat  transfer  coefficient  is  constant  throughout  the  heat  exchanger,  then  we  can 
derive  Eq.  (1)  by  integrating  the  local  differential  heat  flux  equation  from  end  to  end.  Because 
the  transport  properties  of  the  carbon  dioxide  are  temperature  dependent,  the  local  carbon 
dioxide  heat  transfer  coefficient  will  vary  in  the  heat  exchanger,  and  U defined  in  Eq.  (1)  will  not 
always  equal  the  local  value.  The  local  heat  flux  will  also  vary  along  the  tube  length. 

U can  be  further  decoupled  into  the  carbon  dioxide  convection,  tube  conduction,  and  water 
convection  components  through  the  assumption  of  linear  addition  of  the  heat  transfer  modes: 


U 


r r r 1 
• + — ^ln  — + ' 


^C02  ^ss  ri  ro  hf{20 


(2) 


where  hC02  = 
n 

kss 

r0 

hmo  = 


carbon  dioxide  heat  transfer  coefficient; 
inner  tube  inside  radius; 

thermal  conductivity  of  the  stainless  steel  of  the  tube; 
inner  tube  outside  radius;  and 
water  heat  transfer  coefficient. 


The  tube  radius  multipliers  account  for  the  different  heat  transfer  area  for  the  three  terms.  We 
note  that  Eq.  (2)  defines  both  hco2  and  hmo  as  “average”  coefficients  that  are  a single  value  for 
the  entire  tube  length.  In  order  to  determine  hco2  as  a function  of  position  along  the  tube,  we 
would  need  to  know  the  heat  flux  as  a function  of  position,  or  the  specific  enthalpy  of  either  the 
water  or  carbon  dioxide  as  a function  of  position.  The  value  of  hmo  was  measured  in  the 
apparatus  prior  to  performing  the  experiments  of  this  report,  by  operating  the  test  section  with 
hot  water  in  the  inner  tube  and  cold  water  in  the  annular  gap.  hmo  was  a function  of  water  flow 
rate  and  temperature,  which  was  accounted  for  in  its  experimentally  derived  correlation  equation. 
Typical  values  were  7000  to  1 1 000  W/(m2K). 


We  calculated  Q from  the  measurements  of  the  heat  released  by  the  flowing  carbon  dioxide  and 
the  heat  absorbed  by  the  water: 


Q ~ QcOl 


~ mC02  Oo  H )c02  ’ 


(3) 


Q — Qh20 


~ mH20  (h  *0  ) H 20  +Qi 


loss  ’ 


(4) 


where  mco2  = 
Mmo  = 
i — 
Qloss 


carbon  dioxide  mass  flow  rate  (measured); 

water  mass  flow  rate  (measured); 

specific  enthalpy  at  location  L or  0 for  water  or  CO2;  and 

heat  loss  from  water  through  tube  insulation  to  the  room. 


The  specific  enthalpy  was  calculated  from  thermodynamic  equations  of  state  for  the  water  and 
carbon  dioxide  at  the  measured  temperature  and  pressure,  P : 


i = i(T , P) . 


(5) 
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The  equation  of  state  of  Span  and  Wagner  (1996)  was  used  for  the  carbon  dioxide,  and  the  NIST 
Steam  Tables  were  used  for  water  (Gallagher  and  Haar,  1988).  Q!oss  was  measured  as  a function 
of  water  temperature  by  operating  the  water  flow  loop  without  carbon  dioxide  cooling  (Qh20  = 0 
in  Eq.  (4)).  It  ranged  from  0 % to  6 % of  Q,  and  in  some  cases  was  negative  when  the  water  was 
colder  than  ambient. 

To  determine  hco2  at  each  test  point,  the  heat  flow  was  calculated  from  the  measured 
temperatures,  pressures,  and  flow  rates  using  Eqs.  (3),  (4),  and  (5).  QH20  was  used  in  Eq.  (1) 
along  with  the  measured  temperatures  to  calculate  U.  The  value  of  hco2  was  then  calculated 
from  Eq.  (2)  by  subtracting  the  tube  conduction  and  water  convection  terms  from  U. 

In  the  results  that  follow  in  Sec.  4,  we  present  hco2  as  a function  of  the  average  heat  flux  and  the 
carbon  dioxide  mass  flux.  Their  definitions  are: 


Qa  = Q/Aht  = heat  flux; 

(6) 

W = m/Af  = mass  flux. 

(7) 

Af  \s  the  flow  normal  area  for  the  CO2.  We  also  calculated  the  Nusselt  number  (Nu),  along  with 
the  Reynolds  number  (Re)  and  Prandtl  number  (Pr)  for  the  carbon  dioxide  at  each  test  condition. 


„ hD 

(8) 

Nu  = — — , 

k 

Re=?VD, 

(9) 

* 

11 

£ 

(10) 

where  k 

= thermal  conductivity; 

V 

= velocity; 

P 

= density; 

cp 

= specific  heat;  and 

P 

= viscosity. 

Because  the  temperature  of  the  carbon  dioxide  changed  from  the  inlet  to  the  outlet  and  from  the 
tube  center  to  the  wall,  the  fluid  properties  also  changed  throughout  the  flow.  The  fluid 
properties  were  evaluated  at  Tb  (the  average  of  the  CO2  inlet  and  outlet)  in  the  data  analysis.  The 
density  and  specific  heat  were  calculated  from  Span  and  Wagner  (1996),  while  the  viscosity  and 
thermal  conductivity  were  calculated  from  Vesovic  et  al.  (1990).  The  functions  of  Yesovic  et  al. 
include  the  enhancement  of  both  the  thermal  conductivity  and  the  viscosity  around  the  critical 
point. 
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3.3  Experimental  uncertainly 

A summary  of  the  standard  and  combined  standard  uncertainties  is  listed  in  Table  2. 

Uncertainties  in  the  carbon  dioxide  property  functions  and  equation  of  state  were  taken  as  ± 4 % 
at  the  95  % confidence  interval,  based  on  Vesovic  et  al.  (1990)  and  Span  and  Wagner  (1996). 
Combined  standard  uncertainties  in  the  calculated  quantities,  such  as  Q,  U,  hCo2,  Nu,  Re,  and  Pr, 
were  calculated  using  the  “law  of  propagation  of  uncertainty”  as  stated  in  Taylor  and  Kuyatt 
(1994).  This  method  involves  writing  a Taylor  series  approximation  of  those  quantities  and 
applying  the  partial  derivatives  of  the  quantity  with  respect  to  each  variable  times  the  standard 
uncertainty  of  the  variable.  Derivatives  of  the  thermodynamic  and  transport  properties  were 
approximated  with  finite  differences. 

Multiplying  the  values  in  the  Table  by  a coverage  factor  of  2 gives  a 95  % confidence  interval 
that  the  actual  uncertainty  is  less  than  or  equal  to  the  stated  uncertainty.  The  Table  also  lists  the 
major  sources  of  uncertainty  for  each  parameter.  We  can  evaluate  the  quality  of  the 
measurements  by  comparing  the  combined  standard  uncertainty  of  the  heat  balance  error  (Qmcr 
Qcoi)!Qco2  to  the  experimentally  measured  heat  balance  error,  Qerr.  The  experimental  Qerr 
ranged  from  -40  % to  +8.8  %.  Settings  where  the  heat  balance  error  was  larger  than  ±5  % 
corresponded  to  conditions  with  either  the  inlet  or  outlet  CO2  temperature  within  0.2  °C  of  Tm. 
Because  of  the  high  specific  heat  at  Tm,  the  uncertainty  in  the  measured  temperature  produced  a 
large  uncertainty  in  the  enthalpy  and  therefore  the  CO2  heat  flow  at  those  settings.  Twice  the 
calculated  standard  uncertainty  closely  mirrored  the  measured  Qerr  for  those  settings.  If  we 
eliminate  12  of  the  142  settings  with  the  inlet  or  outlet  temperatures  close  to  Tm,  twice  the 
calculated  standard  uncertainty  in  Qen  is  4.1  % to  9.4  %.  For  those  130  settings,  the  measured 
Qerr  ranges  from  —4.7  % to  +4.3  %,  with  a standard  deviation  of  1 .0  %.  Because  the  measured 
uncertainties  in  Qerr  are  less  than  or  equal  to  twice  the  calculated  standard  uncertainties,  we  are 
confident  that  we  have  estimated  the  measured  uncertainties  properly,  and  that  the  uncertainties 
for  the  remaining  calculated  quantities  have  been  estimated  conservatively.  We  have  used  Qmo 
for  Q in  Eq.  (1)  due  to  the  error  produced  by  the  uncertainty  in  the  measured  CO2  temperature  on 
the  CO2  enthalpy. 

Combined  standard  uncertainties  in  hco2  ranged  from  1 .3  % when  the  magnitude  of  hco2  was 
low,  to  9.8  % when  its  magnitude  was  high.  The  largest  source  of  the  uncertainty  in  hco2  was  the 
uncertainty  in  hu20\  at  high  values  of  hco2  the  uncertainty  in  hco2  was  magnified  by  the  ratio  of 
hco^hmo  due  to  the  subtraction  of  the  water-side  heat  transfer  from  the  overall  measured  heat 
transfer  (Eq.  2).  The  combined  standard  uncertainty  in  Numeas  ranged  from  2.4  % to  10.1  %. 

4.  Results  of  experiments 

4. 1 Effects  of  experimental  parameters  on  heat  transfer  coefficient 

Effect  of  W,  Qa , and  P with  Th  > Tm.  The  measured  heat  transfer  coefficient,  hco2,  for  the 
supercritical  carbon  dioxide  from  data  sets  1 to  4 is  plotted  in  Figs.  5 to  7.  For  all  of  these  data 
sets,  the  average  temperature,  Tb,  is  greater  than  Tm.  The  figures  show  representative  data  points 
from  the  experimental  settings,  rather  than  an  average  of  the  multiple  scans  at  each  setting.  Each 
of  these  figures  present  hco2  with  respect  to  mass  flux  (W)  for  different  pressures,  with  heat  flux 
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(Qa)  held  approximately  constant.  The  pressure  is  expressed  as  a ratio  P/Pc.  The  successive 
figures  are  for  increasing  values  of  heat  fluxes.  Two  trends  can  be  seen  on  Fig.  5,  the  lowest  heat 
flux  tested.  First,  hco2  increased  as  the  mass  flux  increased  for  all  pressures.  The  anomalous 
drop  in  hco2,  on  the  P/Pc  =1.05  curve  at  500  kg/(m2s)  mass  flux,  is  due  to  its  higher  Tb  than  the 
other  points  on  the  curve  (we  will  demonstrate  the  effect  of  Tb  in  data  set  6).  The  trend  of  hco2 
increasing  with  W would  also  occur  for  a constant  property  (non-supercritical)  fluid.  Second,  as 
the  pressure  increased,  hco2  decreased  for  all  values  of  the  mass  flux.  At  this  low  QA,  the  heat 
transfer  coefficient  was  highest  at  the  lowest  pressure  (close  to  the  critical)  and  at  high  mass  flux, 
which  also  corresponded  to  high  Reynolds  number. 

Examining  Figs.  6 and  7 for  higher  values  of  the  heat  flux  shows  that  an  increase  in  W always 
increased  hco2 • As  QA  was  increased,  a decrease  in  pressure  produced  a smaller,  yet  still 
significant,  increase  in  hco2.  When  several  sets  of  data  with  the  same  pressure  and  differing  QA 
are  compared,  it  is  seen  that  increasing  the  heat  flux  decreased  hco2  when  the  pressure  was  low, 
but  had  little  effect  on  hco2  when  the  pressure  was  high.  For  a constant  property  fluid,  the  heat 
transfer  coefficient  would  not  change  as  the  heat  flux  varied.  Hence  the  supercritical  fluid 
behaved  more  like  a constant  property  fluid  at  the  high  pressures. 

Effect  of  P/Pr  very  close  to  1 .0.  In  data  set  5,  we  varied  P/Pc  from  1.0034  to  1.125  at 
approximately  constant  flow  rate,  for  four  values  of  Qa.  The  results  for  hco2  are  plotted  in  Fig.  8. 
For  high  QA,  a modest  increase  in  hco2  was  obtained  for  decreasing  the  pressure.  For  the  lowest 
Qa,  hco2  approximately  doubled  as  P/Pc  decreased  from  1.12  to  1.01.  Except  for  the  data  at  P/Pc 
around  1.12,  Tb  > Tm  for  all  the  points  plotted  on  the  figure.  The  slight  decrease  in  hco2  at  the 
lowest  pressure  is  less  than  the  experimental  uncertainty. 

Effect  of  Th  close  to  Tm.  In  data  set  6,  we  varied  Tb  - Tm  from  about  -3  °C  to  +3  °C  at  constant 
pressure  and  approximately  constant  flow  rate.  Fig.  9 shows  the  results  for  four  values  of  Qa- 
For  all  values  of  Qa,  hco2  reached  a peak  when  Tb-Tm  = 0 °C,  although  the  peak  was  much 
larger  for  the  low  heat  flux.  This  figure  also  shows  an  important  relationship  between  the  heat 
flux  and  the  operating  temperature:  for  Tb-Tm>  - 1 °C,  increasing  Qa  caused  hco2  to  decrease, 
whereas  for  Tb  - Tm  < -1  °C,  increasing  Qa  caused  hco2  to  increase.  This  change  in  behavior  with 
temperature  was  not  observed  when  the  carbon  dioxide  was  heated  instead  of  cooled;  for  that 
configuration,  lower  heat  flux  always  improved  the  heat  transfer  coefficient  (Olson,  1999). 

Effect  of  W and  Qa  with  Th  < Tm.  Fig.  10  plots  the  results  of  the  experiments  from  data  set  7 
where  we  varied  W and  Qa  at  P/Pc  = 1 .05  and  Tb-Tm~  -4  °C.  As  we  observed  from  data  set  6, 
hco2  increased  as  Qa  increased.  This  occurred  throughout  the  range  of  flowrates  tested.  hco2 
also  increased  as  W increased. 

4.2  Comparisons  of  the  data  with  predictions  from  the  literature 

The  experimental  results  for  the  Nusselt  number  were  compared  to  several  correlations  for 
turbulent  flow  from  the  literature,  developed  for  either  constant  property  flow  or  supercritical 
flow.  The  first  was  the  constant  property  correlation  of  Petukhov  as  modified  by  Gnielinski 
(1976).  According  to  Rohsenow  et  al.  (1985),  this  is  the  most  accurate  correlation  for  constant 
property  turbulent  flow  in  a tube.  It  has  a stated  accuracy  of  10  % over  2300  < Re  < 5 x 106  and 
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0.5  <Pr<  2000: 


NuCpb  ~ 


f/L  {Re  - 1 OOOlPr  r , , 

/27T-TT73— b(D!Lr\. 


1 + 12.7 


(Pr213  -l) 


(11) 


The  term  in  the  brackets  accounts  for  entry  length  effects.  The  fluid  properties  were  evaluated  at 
the  Tb.  The  friction  factor,/  is  given  by  the  Karman-Nikuradse  correlation: 


= 4.0  log(/te  • y[f )-  0.4 


(12) 


A common  supercritical  correlation  is  the  one  proposed  by  Krasnoshchekov  and  Protopopov 
(1966).  Jackson  and  Hall  (1978)  recommended  this  correlation  for  both  carbon  dioxide  and 
water,  after  they  extensively  reviewed  a number  of  supercritical  correlations  in  the  literature. 

The  correlation  was  meant  to  apply  to  heating  of  supercritical  fluids  when  the  wall  heat  flux  was 
constant.  We  earlier  found  (Olson  and  Allen,  1998)  that  this  correlation  could  be  applied  to  the 
test  section  of  the  present  work  undergoing  convective  heating  boundary  conditions,  and  it  was 
hoped  it  could  also  be  applied  to  convective  cooling  boundary  conditions  to  have  a “universal” 
correlation.  This  correlation  accounts  for  property  gradients  between  the  core  fluid  and  the  wall 
by  applying  a density  ratio  and  a specific  heat  ratio  to  the  constant  property  correlation  (Eq.  1 1): 


Nu-scb  ~ NuCPb 
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The  subscript  “w”  signifies  the  variable  is  evaluated  at  the  conditions  of  the  wall;  the  subscript 
“6”  signifies  an  evaluation  at  the  bulk  flow  conditions  of  the  CO2.  The  wall  temperature  was 
calculated  from  the  linear  addition  of  the  heat  transfer  modes  (Eq.  (2)).  The  exponent  n in  Eq.  13 
depends  on  Tw,  Tb,  and  Tm: 
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For  our  data,  n always  equaled  0.4,  since  TJTm  < 1 .0.  As  noted  by  Pitla  et  al.  (1998),  one  of  the 
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few  correlations  in  the  literature  for  supercritical  cooling  is  that  proposed  by  Krasnoshchekov  et 
al.  (1969),  which  was  developed  from  experimental  cooling  data  on  carbon  dioxide  in  a 
horizontal  tube,  2.22  mm  ID  and  150  mm  long  {L/D  = 67.6).  The  pressure  for  their  experiments 
was  8.0  to  12.0  MPa.  Their  correlation  is: 


Nu  SC  w — NuCPw 


f /O 

p 

1 

1 ^ 

[pb ) 

l CP'W  J 

(15) 


Note  the  similarity  in  form  to  Eq.  13.  Nucpw  is  a constant  property  correlation  of  exactly  the 
same  form  as  Eq.  11,  but  with  all  fluid  properties  evaluated  at  the  wall  temperature.  The 
exponent  p varies  from  0.3  to  0.8,  and  is  a function  of  P/Pc . The  exponent  q is  a function  of  both 
the  specific  heat  ratio  and  P/Pc,  and  varies  from  0.4  to  1 .0.  The  reader  is  referred  to 
Krasnoshchekov  et  al.  (1969)  for  the  p and  q functions. 

Table  3 compares  the  experimental  Nusselt  number,  Numeas,  to  that  predicted  by  the  two  constant 
property  correlations  ( Nucpb , Nucpw ) and  the  two  supercritical  correlations  ( Nuscb , Nuscw)-  For 
each  data  set,  the  extremes  in  the  deviation  between  the  predicted  and  measured  value,  along 
with  the  average  and  standard  deviation  of  the  difference,  are  given  for  each  data  set.  The 
combined  standard  uncertainty  in  the  measurement  is  given  for  comparison.  For  Tb  > Tw,  (data 
sets  1 to  4),  both  constant  property  correlations  under-predicted  the  measured  heat  transfer,  with 
the  wall-based  correlation  performing  worse  than  the  bulk  fluid  correlation.  The  amount  of 
under-prediction  was  more  than  the  combined  standard  uncertainty  of  the  measurement.  The 
wall-based  supercritical  cooling  correlation,  Nuscw , over-predicted  the  measured  values,  also  by 
more  than  the  experimental  uncertainty. 

Remarkably,  the  Nuscb,  developed  for  supercritical  heating , predicted  the  data  for  cooling  when 
Tb  > Tm.  Figure  1 1 shows  the  percent  difference  between  the  correlation  and  measurement 
plotted  against  Tb  - Tm.  The  difference  falls  within  the  10  % uncertainty  band  of  the  measured 
Nu  for  Tb-Tm>  0.  Once  Tb-Tm<  0,  the  differences  become  large.  Analyzing  1200  of  the  1564 
total  points  where  Tb-Tm>  0,  the  average  deviation  between  Nuscb  and  Numeas  was  -1.5%,  and 
the  standard  deviation  was  5.3  % (compared  to  the  range  of  one-sigma  combined  standard 
uncertainty  in  Numeas  of  2.4  % to  10.1%). 

For  the  temperature  range  of  Tb  < Tm  (data  set  7 and  some  settings  from  sets  5 and  6),  the 
correlation  with  the  closest  agreement  to  the  experimental  data  was  the  bulk-based  constant 
property  correlation,  Nucpb ■ It  still  under-predicted  the  data  by  1 1 % to  47  %,  which  was  greater 
than  the  experimental  uncertainty.  All  of  the  other  correlations  also  under-predicted  the  data  for 
Tb  < Tm.  The  Nucpb  correlation  was  unable  to  predict  the  effect  of  the  higher  heat  flux  increasing 
the  heat  transfer  coefficient,  and  also  became  worse  as  the  mass  flow  rate  decreased. 

Figure  12  shows  the  percent  difference  between  Nuscb  and  Numeas  as  a function  of  pressure  ratio, 
for  Tb-Tm>  0.  Although  the  differences  increased  as  P/Pc  approached  1,  they  were  still  within 
two-sigma  of  the  combined  standard  uncertainty  of  the  measurement.  The  correlation  can  be 
confidently  used  over  the  pressure  range  we  tested. 
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5. 


Summary  and  conclusions 


In  this  work  we  have  measured  the  heat  transfer  coefficient  in  turbulent,  supercritical  carbon 
dioxide  flowing  in  a horizontal  tube.  The  carbon  dioxide  was  cooled  by  water  flowing 
countercurrent  in  the  annular  gap  outside  the  tube  containing  the  carbon  dioxide.  This 
convective  boundary  condition  is  the  same  as  in  a heat  exchanger  in  which  the  supercritical 
carbon  dioxide  exchanges  heat  with  another  fluid,  such  as  the  “gas  cooler”  in  a refrigeration 
cycle.  Very  little  experimental  data  exists  in  the  literature  for  cooling  of  a supercritical  fluid, 
particularly  for  tubes  of  diameters  about  1 cm. 

Our  measurements  showed  that  four  parameters  influenced  the  heat  transfer  coefficient  during 
supercritical  cooling:  mass  flow  rate,  cooling  heat  flux,  system  pressure,  and  average  fluid 
temperature,  Tb  (in  relation  to  the  pseudocritical  temperature,  Tm).  The  heat  transfer  coefficient 
always  increased  with  increasing  flow  rate,  as  would  also  be  expected  for  a constant  property 
fluid.  Pressures  close  to  the  critical  produced  high  heat  transfer  coefficients,  with  the  heat 
transfer  coefficient  decreasing  as  pressure  increased.  As  Tb  was  increased  from  below  to  above 
Tm,  the  heat  transfer  coefficient  increased,  reached  a peak,  and  then  decreased  again.  The  effect 
of  heat  flux  on  the  heat  transfer  coefficient  depended  on  Tb  compared  to  Tm\  for  Tb~Tm<- 1 °C, 
an  increase  in  Qa  produced  higher  heat  transfer  coefficients,  while  for  Tb~Tm>  -1  °C,  an 
increase  in  Qa  reduced  the  heat  transfer  coefficient.  For  similar  values  of  these  four 
experimental  parameters,  the  heat  transfer  coefficient  under  cooling  always  exceeded  its  value 
with  heating  (Olson  and  Allen,  1998;  Olson  1999),  at  times  by  a factor  of  two  or  three. 

Several  heat  transfer  correlations  from  the  literature  were  compared  to  the  measured 
experimental  data.  None  were  found  to  adequately  predict  the  measurements  over  the  entire 
range  of  the  experimental  parameters.  A correlation  developed  for  constant  heat  flux  heating 
(Eq.  13)  predicted  the  data  to  within  experimental  uncertainty  when  Tb-Tm>0  °C,  but  under- 
predicted the  data  for  Tb~Tm<0  °C.  This  is  the  same  correlation  that  predicted  all  the  heating 
data  on  this  same  tube.  A correlation  developed  for  supercritical  cooling  (Eq.  15)  over-predicted 
the  data  for  Tb~Tm>  0 °C,  and  under-predicted  the  data  for  Tb~Tm<0  °C.  Two  constant 
property  correlations  under-predicted  the  data  for  nearly  all  conditions,  although  one  of  them 
( Nucpb ) had  the  lowest  errors  of  all  for  Tb~Tm<0  °C. 

All  of  the  correlations  were  developed  for  “local”  heat  transfer  measurements,  whereas  our  heat 
transfer  coefficient  is  defined  as  an  overall  measurement  (using  only  the  inlet  and  outlet 
temperatures  of  the  carbon  dioxide).  It  would  be  useful  to  verify  whether  the  heat  transfer 
equations  could  be  integrated  from  the  heat  exchanger  inlet  to  the  outlet,  with  the  heat  transfer 
coefficient  expressed  by  the  correlations,  to  produce  the  experimentally  measured  conditions. 
This  would  help  determine  if  the  differences  expressed  in  Table  3 are  a result  of  the  differing 
definitions  of  the  heat  transfer  coefficient,  or  truly  represent  the  expected  fundamentally  different 
physics  of  the  supercritical  flow. 
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Table  1.  Summary  of  test  section  parameters  and  experimental  conditions. 
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Boldface  type  indicates  the  parameters  that  were  varied  in  the  data  set. 


Table  2.  Summary  of  standard  uncertainties  for  measurements  and  fluid  properties; 
combined  standard  uncertainties  for  calculated  parameters. 


Parameter 

Major  Source  of  Uncertainty 

Magnitude  of  uncertainty 

CO2  mass  flow  rate 

Instrument  calibration 

0.075  % + 0.34  kg/h 

Water  mass  flow  rate 

Instrument  calibration 

0.05  % + 0.082  kg/h 

Fluid  temperature 

Instrument  calibration 

0.025  °C 

CO2  absolute  pressure 

Instrument  calibration 

2070  Pa 

CO2  differential  pressure 

Instrument  calibration 

43  Pa 

Water  absolute  pressure 

Instrument  calibration 

2600  Pa 

Water  differential  pressure 

Instrument  calibration 

430  Pa 

CO2  specific  enthalpy  and 
specific  heat 

Equation  of  state 

2% 

CO2  viscosity  and  thermal 
conductivity 

Function  accuracy 

2% 

Heat  balance  error,  all  settings 

Temperature  measurement, 

CO2  specific  enthalpy  function 

2.1  % to  19.7% 

Heat  balance  error,  exclude  12 
settings  where  inlet  and  outlet 
CO2  temps  within  0.2  °C  of  Tm 

Temperature  measurement, 

CO2  specific  enthalpy  function 

2.1  % to  4.7% 

Water  heat  transfer  coefficient 

Curve  fit  to  NIST  calibration 
data 

4.0%  to  5.1  % 

CO2  heat  transfer  coefficient 

Water  heat  transfer  coefficient 

1.3%  to  9.8% 

CO2  Nusselt  number 

CO2  heat  transfer  coefficient, 
thermal  conductivity 

2.4%  to  10.1  % 

CO2  Reynolds  number 

CO2  viscosity 

2.0  % to  2.5  % 

CO2  Prandtl  number 

Temperature  measurement, 

CO2  specific  heat 

3.5  % to  7.7  % 

Multiplying  the  Table  values  by  a coverage  factor  of  2 will  yield  a 95  % confidence  interval. 
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Table  3.  Comparison  of  the  Nusselt  number  predicted  by  four  correlations  to  measured  values. 
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Figure  1.  Properties  of  carbon  dioxide  on  the  7.8  MPa  isobar  from  5 °C  to  55  °C. 
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Figure  2. 


Carbon  dioxide  flow  loop  of  the  NIST  supercritical  heat  transfer  facility. 


Cold 


Figure  3. 


Water  flow  loop  of  the  NIST  supercritical  heat  transfer  facility. 
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Figure  4.  Counterflow  heat  exchanger  test  section. 
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Figure  5.  CO2  heat  transfer  coefficient  (hcoi)  as  a function  of  mass  flux  (W)  for 
several  values  of  pressure  (P/Pc)  at  a heat  flux  (QA)  of  20  kW/m2  to  22 
kW/m2.  Tb  > Tm. 
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Figure  6.  CO2  heat  transfer  coefficient  {hcoi)  as  a function  of  mass  flux  (W)  for 
several  values  of  pressure  (P/Pc)  at  a heat  flux  (QA)  of  42  kW/m2  to  44 
kW/m2.  Tb  > Tm. 
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Figure  7.  CO2  heat  transfer  coefficient  ( hco2 ) as  a function  of  mass  flux  (W)  for 

several  values  of  pressure  (P/Pc)  at  a heat  flux  (Qa)  of  61  kW/m2  to  65 
kW/m2.  Tb  > Tm. 
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Figure  8.  CO2  heat  transfer  coefficient  (hcoi)  as  a function  of  pressure  (P/Pc)  for 
several  values  of  heat  flux  ( Qa ) at  mass  flux  (W)  of  569  kg/(m2s)  to  613 
kg/(m2s). 
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Figure  9.  CO2  heat  transfer  coefficient  (/2C02)  as  a function  of  bulk-to-pseudocritical 
temperature  difference  (7*  - Tm)  for  several  values  of  heat  flux  (Qa). 
Pressure  (P/Pc)  of  1 .05  and  mass  flux  (W)  of  540  kg/(m2s). 
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Figure  10.  CO2  heat  transfer  coefficient  (hcoi)  as  a function  of  mass  flux  (W)  for 

several  values  of  heat  fltix  (Qa)-  Pressure  (P/Pc)  of  1.05  and  Tb  < Tm. 
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SCb  meas 


Figure  1 1 . Difference  in  Nusseit  number  between  Nu$cb  correlation  and  Numeas  as  a 
' 1 function  of  bulk-to-pseudocritical  temperature  difference  (Tb  - Tm). 
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Figure  12. 


. ~ isV-  ijss&u^.  rri  rormaft  ’i  :•  1 1 

Difference  in  and  Numeas  as  a 

function  of  pressure  (P/Pc)  for  Tb  > Tm. 
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